Exact singlet bond ground states for electronic models 
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We have proposed several ID and 2D electronic models with the exact ground state. The ground 
state wave function of these models is represented in terms of "singlet bond" functions consisting of 
homopolar and ionic configurations. The Hamiltonians of the models include the correlated hopping 
of electrons, pair hopping terms and the spin interactions. One of the ID models demonstrates 
spontaneous breaking of the translational symmetry. We studied also the model describing the 
transition point between the phases with and without an off'-diagonal long-range order. 



The study of strongly correlated electron systems has 
been an important subject in theoretical condensed mat- 
ter physics. In general, the Hamiltonians of these systems 
include many types of interactions, and they are difficult 
to solve. The integrable models provided us a very good 
understanding of the correlation effects in many-body 
systems. Unfortunately, the construction of such models 
is difficult due to the strict conditions for the integrabil- 
ity. In recent years there has been increasing interest in 
studying models where at least the ground state can be 
found exactly The most popular methods for the 

construction of exact ground state are so-called optimal 
ground state (OGS) approach ||^ and the matrix-product 
(MP) method The ground state wave function in 

the MP method is represented by Trace of a product of 
matrices describing single-site states. This ground state 
is optimal in the sense that it is the ground state of each 
local interaction. This method allows to construct a large 
class of spin models. The similar spirit has been used in 
the OGS method for construction of the electronic mod- 
els with the special ground states. 

In this Letter we propose new ID and 2D models of in- 
teracting electrons with the exact ground state. We note 
that our models have ground states which are very dif- 
ferent from those constructed in the OGS approach. The 
ground state wave function of our models is expressed 
in terms of the two-particle "singlet bond" (SB) function 
located on sites i and j of the lattice: 



)|0), (1) 



where cf^, ci^a are the Fermi operators and x is an ar- 
bitrary coefficient. The SB function is the generalization 
of the RVB function Q including ionic states. The pres- 
ence of the ionic states is very important from the phys- 
ical point of view because, as a rule, the bond functions 
contain definite amount of the ionic states as well. 

It is known a set of ID and 2D quantum spin models 
the exact ground state of which can be represented in 
the RVB form |^-^ . It is natural to try to find electronic 
models with exact ground state at half-filling formed by 
SB functions in the same manner as for above mentioned 
spin models. The electronic models of these types include 
the correlated hopping of electrons as well as the spin 
interactions and pair hopping terms. 



The model with dimerization.- As the first example 
we consider the ID electronic model with the two-fold 
degenerate ground state in the form of the simple product 
of SB dimers, similarly to the ground state of the well- 
known spin-i Majumdar-Ghosh model 1^. For the half- 
filling case the proposed ground state wave functions are: 



and 



^-1 [l,2][3,4]...[iV-l,iV] 



= [2,3][4,5]...[7V, 1] 



(2) 



(3) 



In order to find the Hamiltonian for which the wave 
functions (^ and (^ are the exact ground state wave 
functions, we represent the Hamiltonian as a sum of lo- 
cal Hamiltonians hi defined on three neighboring sites 
(periodic boundary conditions are supposed): 



N 



H 



(4) 



The basis of three-site local Hamiltonians hi consists 
of 64 states, while only eight of them are present in 'i'l 
and \l/2. These 8 states are 



[i,i + l]ipi+2, (pi[i + l,i + 2], 



(5) 



where ipi is one of four possible electronic states of z-th 
site: |0),, |t)„ U)„ |2),. 

The local Hamiltonian hi for which all the functions (j^) 
are the exact ground state wave functions can be written 
as the sum of the projectors onto other 56 states \xk) 



hi = ^Xk \Xk) iXk 



(6) 



where are arbitrary positive coefficients. This means 
that the wave functions ^'i and ^'2 are the ground states 
of each local Hamiltonian with zero energy. Hence, ^1 
and ^'2 are the ground state wave functions of H with 
zero energy. 

In general case, the local Hamiltonian hi is many- 
parametrical and depends on parameters Afc and x. We 
consider one of the simplest forms of hi including the cor- 
related hopping of electrons of different types and spin 
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FIG. 1. The lattice of the Shastry - Sutherland model 

interactions between nearest- and next-nearest neighbor 
sites: 

hi = 2 — X (ti^i+i + <i+i,j+2) 

+ (a;2_(l + a;2)(l-n,+i)2)T,^,+2 
1 — 

+ 8 (Si ■ Si+i + Si+i • S,;4-2 + Si • Si+2), (7) 

where 

a 

(7 

and Si is the SU{2) spin operator. 

Each local Hamiltonian hi is a non-negatively defined 
operator at |a;| < 1. The following statements related to 
the Hamiltonian (|^) are valid: 

1. The functions (|) and (||) are the only two ground 
state wave functions of the Hamiltonian at A^e = ^ 
{Ne is the total number of electrons) . They are not or- 
thogonal, but their overlap is ^ at iV 3> 1. 

2. The ground state energy Eo{Ne/N) is a symmetrical 
fimction with respect to the point N^/N = 1 and has a 
global minimum i?o = at A^e/.^ = 1- 

3. The translational symmetry of (0) is spontaneously 
broken in the ground state leading to the dimerization: 

{\'ti,i+l ~ ti+l,i+2\) = 2 

The excited states of the model can not be calculated 
exactly but we expect that there has to be a gap, because 
the ground state is formed by the ultrashort-range SB 
functions. If it is the case, the function Eo{Ne/N) has a 
cusp a.t Ne/N = 1. 

Actually, this model is the fermion version of the Ma- 
jumdar - Ghosh spin model. Moreover, it reduces to the 
Majumdar - Ghosh model at x = and in the subspace 
with Hi — 1. 

For a; = 1 the Hamiltonian (|^) simplifies and takes the 
form: 

H ^ -2^(t,,,+i - 1) - 5]e-^+^T„-+2 (8) 
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FIG. 2. The two-leg ladder model. 



The 2D model- We can easily construct the 2D electronic 
model with the exact ground state which is analogous to 
the Shastry - Sutherland model (Fig.l). The Hamil- 
tonian of this model is: 

Ki+Kk + hlf^, (9) 

where the sum is over all triangles {i, j, fc}, one of which is 
shown on Fig.l. So, each diagonal line belongs to the two 
different triangles. The local Hamiltonians /i^^, acting on 
the diagonal of the triangle fc}, and /lij, hi^k have 
the form (for the sake of simplicity we put a; = 1) 

hik = "2 t,,fe + 4 

hi^k — ii,k ^ ^ Ti k 

It is easy to check that 

hj,k \fi = [K^j + h^^k) Wi [j.k]) = 

All other states of the local Hamiltonian hij+hi^k+hj f. 
have higher energies. Therefore, the ground state wave 
function in the half-filling case is the product of the SB 
functions located on the diagonals shown by dashed lines 
on Fig.l. This model has non-degenerate singlet ground 
state with ultrashort-range correlations. 

The ladder model- Let us now consider electronic 
models with a more complicated ground state including 
different configurations of short-range SB functions. The 
form of these ground states is similar to that for spin 
models proposed in ^ and generalized in ||l^. In the ID 
case our model describes two- leg ladder model (Fig. 2). 
Its ground state is a superposition of the SB functions 
where each pair of nearest neighbor rungs of the ladder 
is connected by one SB. One of possible configurations of 
singlet bonds is shown on Fig. 2. 

The wave function of this ground state can be written 
as: 

= V^^il)g^,r'{2)gp. ■ . . r^{N)grx (10) 

The wave function of this type for spin models has been 
proposed in The functions describes i-th rung 

of the ladder 

V'^^(l) - Ci^2r-lV2^ + Cw\^^2^-l (H) 

with 
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It is easy to see that 
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Therefore, the function ^'^ is a singlet wave function 
depending on two parameters x and cxjci. Actually, this 
form of \l/s is equivalent to the MP form with 4x4 
matrices A\y[i) = gx^^ip^'^ (i). Moreover, at a; = and 
Ci/c2 = — 1 the function reduces to the wave function 
of the well-known AKLT spin-1 model. 

In order to find the Hamiltonian for which the wave 
function ( [lo| ) is the exact ground state wave function, it 
is necessary to consider what states are present on the 
two nearest rungs in the '^s- It turns out that there are 
only 16 states from the total 256 ones in the product 
il)^^{i)g^^ip'^P{i + 1). The local Hamiltonian hi acting on 
two nearest rungs i and \ can be written in the form of 
with the projectors onto the 240 missing states. The 
total Hamiltonian is the sum of local ones (^. The ex- 
plicit form of this Hamiltonian is very cumbersome and, 
therefore, it is not given here. 

The correlation functions in the ground state (|l^) can 
be calculated exactly in the same manner as it was done 
for spin models IJnj. It can be shown that all of correla- 
tions exponentially decay in the ground state. We expect 
also that this model has a gap. 

This method of construction of the exact ground state 
can be generalized also to 2D and 3D lattices 

©• Fol- 

lowing pOl , one can rigorously prove that the ground 
state of these models is always a non-degenerate singlet. 

ID models with the giant spiral order- There is one 
more spin-^ model with an exact ground state of the 
RVB type H. Its Hamiltonian has the form 



H 



'i+l 



'i+2 



(12) 



This model describes the ferromagnet-antiferromagnet 
transition point ||ll[ . The exact singlet ground state can 
be expressed by the combinations of the RVB functions 
(i, j) distributed uniformly over the lattice points: 



(13) 



where the summation is done over all combinations of 
sites under the condition that i < j, k < I, m < n . . .. 
The spin correlations in the singlet ground state show 
giant spiral structure 

The analog of the wave function (|l^) in the SB terms 

is: 



i<j.. 



(-1)^ [t,j][k,l][m,n] 



(14) 



where P = {i,j, k,l, . . .) is the permutation of numbers 
(1,2,... N). It is interesting to note that the singlet wave 
function (|lj) can be also written in the MP form but with 
an infinite size matrices fl^ . 

In order to find the Hamiltonian for which the wave 
function (|l^) is the exact ground state wave function, let 
us consider what states are present on the two nearest 
sites in the ^o- It turns out |13[ t hat there are only 9 
states from the total 16 ones in (|14D. They are 



ITT), lii), ITi + iT), 120 - 02} 

Iti- it) -a; 120 + 02), IT 0-0 T), 
IT 2-2 T), liO-OT), IT 2-2 i) 



(15) 



The local Hamiltonian /li.i+i can be written as the sum 
of the projectors onto the 7 missing states (^. At \x\ > 1 
the most simple form of this total Hamiltonian is: 



H 



N 

E 



T, 



iA+l H 

X 



4 S, • S 



4 x'^ 
—Vt ■ m+i + 4 — 



We use here 77 operators: 



(16) 



which form another SU{2) algebra |12 ij], and ryi ■ 772 
is a scalar product of pseudo-spins r]i and 772- We note 
that this Hamiltonian commutes with S^, but does not 
commute with 77^. It can be prove d [p^ that only three 
multiplets are the ground states of the singlet state 
(|T^), the trivial ferromagnetic state S — N/2 and the 
state with S = N/2 - 1. 

The norm and the correlators of the electronic model 
(M) in the singlet ground state are exactly calculated 
p3[. For example, the norm of (n4) is: 



/T IT X d'^ / l + cosh(a;0 



The correlators at ^ 1 are 



{eta = 0(^ 



I 'li+ll 



(S,S,+,)- Jcosl"^) (17) 



The correlator (vtVi+i) which determines the off- 
diagonal long-range order (ODLRO) [|l2| vanishes in the 
thermodynamic limit. At the same time the spin-spin 
correlations have a spiral form, and the period of the 
spiral equals to the system size as in the spin model ([l^ ) . 
Another electronic model can be obtained by making 



the canonical transformation cf 



4a ^"^^ %i 
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As a result of this transformation, the SB function (|l|) 
becomes: 

and the wave function ( p^ changes to 

*o- (19) 

i<j... 



The function ( |19D for > 1 is the exact ground state 
wave function of the transformed Hamiltonian: 



N 



1=1 

^S, • S.+i + 4 ^—-^StSt+, + 1 ) , (20) 



where 



x^ 



Similarly to the original spin model (12) [ p]JT^ the last 
two electronic models ([l^),(po|) describe the transition 
points on the phase diagram between the phases with and 
without a long-range order [ferromagnetic for the model 
( |T6| ) and off-diagonal for the model (|2^)]. Therefore, we 
suggest the formation of the ground state with giant spi- 
ral order (ferromagnetic or off-diagonal) as a probable 
scenario of the subsequent destruction of the ferromag- 
netism and superconductivity. 

In summary, we have constructed electronic models 
with an exact ground state. The ground state wave 
function of these models is built from SB functions in 
the same manner as the well-known RVB ground states 
of spin models. We have considered three types of SB 
ground states. However, we note that the proposed ap- 
proach can be generalized for the construction of other 
models with the ground states of more complicated SB 
forms. 

Authors are grateful to Max-Planck-Institut fur Physik 
Komplexer Systeme for a kind hospitality. This work was 
supported by RFFR. 



This Hamiltonian commutes with rf' but does not com- 
mute with S^. Therefore, the eigenfunctions of the 
Hamiltonian ( pO| ) can be described by quantum numbers 
•q and 77^ . For the cyclic model ( pO| ) the states with three 
different values of 77 have zero energy [as it was for 
the model (p^)]. They include one state with 77 = (19), 
all states with r/ = A^/2: 



|0) 



(21) 



and the states with rj = A^/2 — 1. Therefore, for the case 
of one electron per site (77^ — 0), the ground state of the 
model (20) is three-fold degenerate. 

The correlation functions in the ground states with 
77 = N/2 and 77 = N /2 — \ for the half- filling case coincide 
with each other and at 3> 1 they are: 



N 



The existence of the ODLRO immediately follows from 
the latter equations. The correlation functions in the 
ground state ( p^ ) have similar forms as in Eqs.(p^: 



(23) 



The giant spiral ordering in the last equation implies the 
existence of the ODLRO and, therefore, the supercon- 
ductivity [|l^ in the ground state (p^). 
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